Response function and elementary excitations of a Bose-Einstein condensed gas with 
equal Rashba and Dresselhaus spin-orbit coupling 
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By calculating the dynamic density response function we identify the excitation spectrum of a 
Bose-Einstein condensate with equal Rashba and Dresselhaus spin-orbit coupling at zero temper- 
ature. The propagation of phonons is strongly affected by the coupling and the velocity of sound 
vanishes when one approaches the second-order phase transition between the spin polarized and 
the zero-momentum quantum phases. We also point out the emergence of a roton minimum in the 
excitation spectrum for small values of the Raman coupling, when one approaches the transition to 
the stripe phase. 
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Synthetic gauge fields are a developing field of research 
in atomic physics. They have been the object of recent 
relevant experimental [IHl] and theoretical work fj- 15|. 
giving rise to the occurrence of new quantum phases ex- 
hibiting unique magnetic features. The elementary exci- 
tations of such systems are also expected to exhibit novel 
features [lil - Elj ]. Some of these features have been al- 
ready the object of experimental measurements In 
particular the experiment of Q has shown that the cen- 
ter of mass oscillation of a harmonically trapped Bose- 
Einstein condensate can be deeply affected by the cou- 
pling with the spin degree of freedom, in agreement with 
the predictions of theory (2fjj . 

The purpose of the present work is to study the el- 
ementary excitations and the corresponding behavior of 
the dynamic structure factor of a spin-orbit coupled Bose- 
Einstein condensate (BEC) at zero temperature, by di- 
rect investigation of the response of the gas to a time de- 
pendent perturbation. We explore both the the phonon 
regime of long wavelengths and the region at higher mo- 
mentum transfer where the spectrum exhibits novel fea- 
tures including the occurrence of a roton minimum. In 
ultracold gases the excitation spectrum can be measured 
via two-photon Bragg spectroscopy [22], so our predic- 
tions can be relevant for future experiments on spin-orbit 
coupled BEC's. 

We consider a spin 1 /2 interacting Bose gas of N parti- 
cles enclosed in a volume V, characterized by the Hamil- 
tonian (we set h = m = 1) 



H = h o{i) + H\J d3r ^P n a (r) n/3 {r) (1) 

i a, j5 

where i = 1, • • • , N is the particle index, while a and 
f3 are the spin indexes (t, 4 = i) characterizing the two 
spin states. The single particle Hamiltonian /i is defined 
by 



and is characterized by equal contributions of Rashba [23] 
and Dresselhaus [24[ spin-orbit couplings. f2 is the Ra- 
man coupling accounting for the transition between the 
two spin states, fco the momentum transfer of the two Ra- 
man lasers and 5 is an effective magnetic field providing 
the energy difference between the two single-particle spin 
states. The spin up and down density operators entering 
Eq.© are defined by n±(r) = (1/2) £\ (1 ± cr M ) 5(r-n) 
while g a @ = 4ira a p are the relevant coupling constants in 
the different spin channels, with a a /3 the corresponding s- 
wave scattering lengths. Finally a^, with k — x, y, z, are 
the usual 2x2 Pauli matrices. The Hamiltonian (jXJ[2l) has 
been already implemented experimentally 0, [E| and has 
been recently employed to describe a variety of non triv- 
ial quantum phases in Bose-Einstein condensates 11 , 12| ■ 
It has the peculiar property of violating both parity and 
time-reversal symmetry. 

In the presence of a spin symmetric Hamiltonian (g^ = 
9ii = 9 an d S — 0), the mean field ap pro ach predicts the 
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occurrence of three quantum phases [llj, [12[ , depending 
on the value of the relevant parameters ko, the density 
n, and the interaction parameters 

Gi = n (g + g n ) /4 ; G 2 = n {g - g n ) /4 (3) 

These are the stripe phase (also called phase I) where 
Bose-Einstein condensation (BEC) takes place in a com- 
bination of two plane waves with opposite momenta, the 
single plane wave or spin polarized phase II, where BEC 
takes place in a single plane wave state with momen- 
tum different from zero, and the zero momentum phase 
III where the condensate has zero momentum. The crit- 
ical value f^ 1-11 ) of the Raman coupling characterizing 
the first order transition between the stripe and the 
spin polarized phases can be calculated analytically in 
the small coupling limit G\, Gi <C k\, where one finds 
ftC" 1 ) = 2^^/27/(1 + 27) with 7 = G2/G1. The tran- 
sition between the phases II and III takes instead place 
at the higher value f^ 11 " 111 ) = 2 (fcg - 2G 2 ) , provided the 
condition k% > 4G 2 (1 + 7) is satisfied, and has a sec- 
ond order nature, being characterized by the divergent 
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behavior of the spin polarizability [20j. If G 2 < only 
the phases II and III are available, the stripe phase being 
always energetically unfavorable. 

In order to calculate the dynamic response of the 
system we add the time dependent perturbation V\ = 
_^ e t(qx-ut) _|_ jj c -(- |- ne s i n gi e particle Hamiltonian 

@. The perturbation provides momentum transfer q 
along the most relevant x-direction. The density response 
function is then calculated through the usual definition 
x(q,w) = lim\^ f oSn/(Xe~ lu ' t ) 7 where 5n are the density 
fluctuations induced by the external perturbation. In the 
following we will calculate x(<3S <*0 by solving the time de- 
pendent Gross-Pitaevskii equations. We will restrict the 
analysis to the phases II and III, where the ground state 
density is uniform and the wave function of the conden- 
sate can be written in the simple form 



cos 9 
- sin# 
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In the phase II the ground state wave function (first term) 
is characterized by the condition 9 = a.rccos(k 1 /k )/2 
and the value k x = k oy /l - « 2 /[ 4 ( fc o = 2G 2) 2 ] for the 
momentum where BEC takes place, corresponding to the 
value (rif — ni)/n = fci/fco for the spin polarization. In 
the phase III one instead has fci =0, 9 = 7r/4 and — 
H- 

By calculating the Bogoliubov amplitudes u and v in 
the linear (small A) limit we find the result (near the 
poles one should replace u with uj + iO) 



N [-g 2 m 2 + ai(g)u; + 00(g)] 
w 4 + b 3 (q)uj 3 + b 2 (q)oj 2 + &i(g)w + b (q) 



(5) 



where the values of the various parameters depend on 
whether one is in Phase II or III. In phase II we find 

ao = g 6 /4 - [(3 + 2G 2 ) k\ - 1 + 2G 2 ] g 4 

+ 4(l-2G 2 ) 2 fc 2 g 2 
ai = 4fcig 3 

6 = q 8 /16~ (k 2 1 -G 1 + G 2 )q 6 /2 

+ [kf-2 (3d + 4G 2 + 2GiG 2 + 2G^) k\ 
+ 2G!(l-2G 2 )]g 4 

+ 4(1- 2G 2 ) [G 2 k 2 + G x (1 - 2G 2 ) - 2G 2 ] k\q 2 
b x = k iq 5 - 4 (k 2 - 1 - 2Gi - 2G 2 ) fcig 3 

- 16G 2 (1 - 2G 2 ) (k\ - 1) k iq 
b 2 = - q 4 /2 + 2 (3k 2 -2-G1 + G2) q 2 

+ 4(l-2G 2 )(2G 2 fc 2 -l) 
63 = -4fcig 



In Phase III we instead obtain the results 

a = g 6 /4 + (0 - 1 + 2G 2 ) g 4 + (1 [SI - 2 (1 - 2G 2 )] g 2 
b = g 8 /16 + [fi - 2 (1 - Gi - G 2 )] g 6 /4 
+ [fi 2 -4(l-2Gi-G 2 )fi 

+ 4(l-2G!)(l-2G 2 )]g 4 /4 

+ 2GiO[17-2(l-2G 2 )]g 2 
b 2 = -g 4 /2 - [fi + 2 (1 + Gi + G 2 )] g 2 - fl (fl + 4G 2 ) 
ai = b% = 63 = 

where, to simplify the notation, we have expressed the 
quantities g and k\ in units of k , and the quantities 
f2, Gi and G 2 in units of k 2 . The above equations 
include all the relevant information relative to the fre- 
quency of the elementary excitations (poles of x) , as weu 
as to the dynamic structure factor given, at T — 0, by 
S(q,u) = 7r _1 Imx(g,w) for w > and S(q,uj) — for 
negative to. They also permit to calculate the static 
response function x(q) = x(9: tJ — 0)/-^ yielding, for 
q — > 0, the compressibility x 01 the gas. It is worth 
noticing that the odd terms in uj entering the response 
function identically vanish in the spin unpolarized phase 
III, but survive in the spin polarized phase reflecting 
the lack of parity and time reversal symmetry of the 
ground state wave function in phase II. The condition 
Imx(g,o;) = — Inr\;(— g, — uS) characterizing the imaginary 
part of the response function is always satisfied, but the 
symmetry relationship Imx(g,w) = Imx(— q, w) is not 
ensured in the phase II. Despite this asymmetry and the 
fact that the compressibility exhibits a jump at the tran- 
sition, the frequencies of the elementary excitations are 
always continuous functions of ft. First results for the 
excitation spectrum of the Hamiltonian ([I])-® for small 
and large values of fi, far from the transition between the 
plane wave and single minimum phases, have been re- 
cently discussed in 21| using a hydrodynamic approach. 

The low frequency excitations at small q (sound waves) 
can be easily obtained by solving the equation b 2 tu 2 + 
b\u> + 60 = 0. By setting w = eg we find, in the phase III, 
the result 
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(6) 



which explicitly shows the strong reduction of the sound 
velocity caused by the spin-orbit coupling when one ap- 
proaches the transition to the spin polarized phase II. 
At the transition the sound velocity vanishes and the el- 
ementary excitations exhibit a different q 2 dependence 
along the x-direction. 

In the plane wave, spin polarized phase II the behavior 
of the sound velocity exhibits a further interesting feature 



3 



0.7 
0.6 
0.5 
£ 0.4 
o 0.3 
0.2 
0.1 














- 










c iN\ 











1.5 



2.5 



f2/kn 



FIG. 1: (Color online) Sound velocity as a function of the Ra- 
man coupling for the following choice of parameters: Gi/fco = 
0.2, G2/feo = 0.05. The two sound velocities in the Phase II 
correspond to phonons propagating in the direction parallel 
(c^) and antiparallel (cj^) to k\. The horizontal dashed line 
corresponds to the value \/2Gi = 0.63 ko of the sound veloc- 
ity in the absence of spin-orbit and Raman coupling. The 
vertical dashed lines indicate the Raman frequencies at which 
the I-II and II-III phase transitions take place. 



caused by the lack of symmetry between q and — q 
± _ W 2 (*o - 2G 2 ) [(fc 2 _ 2d - 2G 2 )G 2 fc 2 + fc^j 



fc* - 2G 2 fc 2 



± 



2G 2 fc! (fc 2 - k\) 
2G 2 k\ 



fc 4 



(7) 



where the sign +(— ) indicates sound waves propagating 
parallel (anti-parallel) to k\ . Notice that the asymmetry 
effect in the sound velocity is caused by the spin inter- 
action term G 2 . Setting G 2 = one finds c^J = c^j = 
\j2Gi[l - £l 2 /(4fc;j)]. Also in the phase II the sound ve- 
locity vanishes when one approaches the transition to 
the phase III where fci = 0. The results for the sound 
velocity are shown in FigQ] for a configuration with rel- 
atively large G 2 emphasizing the difference between c^f 
and Notice that the sound velocity, in the absence of 
spin-orbit and Raman coupling, would correspond to the 
value c = \J2G\. This value is asymptotically reached 
only for very large values of Q. 

The strong quenching effect exhibited by the sound 
velocity near the II-III phase transition is a remark- 
able consequence of spin-orbit coupling. This effect 
is surprising especially in the zero momentum phase 
III where the compressibility of the gas is unaffected 
by the coupling and given by x = l/(2Gi). The 
quenching of the sound velocity is the consequence of 
the strong coupling between density and spin density 
in the propagation of sound. It also reflects the fact 
that, as can be inferred from the analysis of Eq.©, 
the phonon modes exhaust the compressibility sum rule 
f^ 00 du)S(q,ui)/(jj but, differently from ordinary super- 
fluids, they give a vanishingly small contribution to the 
/-sum rule f + ™ duiS(a, lo)uj — Nq 2 /2 when one ap- 
proaches the transition [25]. This behavior shares impor- 
tant analogies with the quenching of the frequency of the 



center of mass oscillation, recently explored both experi- 
mentally Q and theoretically [2(|. Also the static struc- 
ture factor S(q) — / °° duiS(q,co) is strongly quenched 
compared to usual BEC's. This results in an enhance- 
ment of the quantum fluctuations of the phase as pre- 
dicted by the uncertainty principle inequality [26| . 

The peculiar behavior of the excitation spectrum in 
the phonon regime can be usefully described also using 
the hydrodynamic formalism where one writes the spin- 
up and spin-down components of the order parameter in 
terms of their modulus and phase 21] . In this case one 
finds four coupled equations instead of two equations as 
in usual BEC's. In the phonon regime of large wave- 
lengths and small frequencies one can safely neglect the 
quantum pressure terms. Furthermore, one finds that 
the phase difference between the two spin components is 
blocked (<p^ — tp^) [27]. By requiring the vanishing of 
the difference d t ipf — d t fi between the equations for the 
two phases, one then derives the non trivial relationship 
(<P = <Pt = <Pl) 



k V x 5ip - k Z 
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between the phase gradient, the density and the spin fluc- 
tuations. In the above equation s — nki/ko is the spin 
density relative to the equilibrium configuration while 
Z = fl/{2k 2 {l - fc 2 /fc 2 ) 3 / 2 ). Eq.© permits to reduce 
the relevant hydrodynamic equations, i.e. the equation 
of continuity 



d t Sn + V ■ (riVStp) - fc V 2 



s Sn 5s 



(9) 



and the equation for the phase 

Si S S Tl Ss 

dtSip + k x \7 x 5(p - k 2 Z~ [-- h — 



2 Gl ^=0, 
n 

(10) 

to a closed set of coupled equations for the density and 
the phase. It is worth pointing out the crucial changes 
caused by the spin orbit term in the equation of continu- 
ity ([3]). These changes reflect the fact that the physical 
current of the system is not simply given by the momen- 
tum operator, but is also affected by the spin variable. 
One can identify the current density operator as j(r) = 
p(r) - k a z (r)i, where p(r) = £\ [p* 6(r - i-;) + H.c] /2 
and cr z (r) = a z,i 5(r— r,) are the momentum and spin 
density, respectively, and i is the unit vector along x. 
The new current operator satisfies the continuity equa- 
tion [H,n(r)} = iV ■ j(r), with n(r) = J2i^( r — r tnc 
density operator. In deriving the equation for the phase 
we have also identified the expression 



(i = 2Gi + 2G 2 



■ (k 2 - 2G 2 Y 



(11) 



for the chemical potential in phase II whose density de- 
pendence is affected by spin-orbit coupling only if G 2 =/= 0. 



4 




In phase III one has instead fi = 2G\ + (k 2 — f2)/2, so 
that fi is continuous at the transition. 

The solutions of the hydrodynamic equations repro- 
duce exactly the results © for the sound velocity along 
the x axis. Concerning the sound velocity along the other 
directions y and z we find that in phase III its value is not 
modified with respect to the value \J2G\ of usual BEC's, 
while it is affected in phase II if G2 7^ 0. 

The hydrodynamic equations further permit to calcu- 
late the relative amplitudes of the density, spin density 
and momentum density oscillations characterizing the 
propagation of sound. If G% = we find the result 



(12) 



where \ = l/(2Gi) is the compressibility and xm is the 
magnetic polarizability, defined by x~m = Nd 2 E / dM 2 = 
k 2 (Z - 1) |!3|, with E and M = JV t - the total energy 
and magnetization of the system, respectively. Equation 
(fT2)) shows that, near the transition between phases II 
and III, where \m diverges, the density fluctuations are 
suppressed and the system behaves like an incompressible 
fluid, despite the fact that the compressibility x is finite. 

A major usefulness of the hydrodynamic equations is 
that they can be easily extended to trapped non-uniform 
configurations. In the simplest G2 = case, correspond- 
ing to G\ — ng/2, where the wave vector k\ is density 
independent, the 3D hydrodynamic equations can be re- 
duced to the compact form: 

d?5n = g [(1 - l/Z) V* (nV x Sn) + V_l (nV ± Sn)} (13) 

with n the Thomas-Fermi density profile given, in the 
presence of harmonic trapping Vho(r) = (oj 2 x 2 + w 2 y 2 + 
lo 2 z 2 )/2, by an inverted parabola in both phases II and 
III: n(r) = [/i — Mio(r)] /g, with /x fixed by the nor- 
malization condition. One can easily check that all the 
solutions holding for usual BEC's [28( still hold in the 
presence of spin-orbit coupling, with the simple replace- 
ment of the trapping frequency w x with w x ^l — l/Z. 
This reproduces the result derived in [2(J for the dipole 
oscillation using a sum rule approach and also shows that 
the frequency of the other hydrodynamic modes involving 
a motion of the gas along the x axis will be consequently 
quenched. 

In the last part of this work we discuss the behavior 
of the excitation spectrum at finite momentum transfer 
q. The poles of Eq.© actually provide two separated 
branches, the lower one approaching the phonon disper- 
sion at small q. A very peculiar feature of the lower 
branch is exhibited in the plane wave spin polarized phase 
II for negative values of q, resulting in the emergence of a 
roton minimum 2l| which becomes more and more pro- 




nounced as one approaches the phase transition to the 
stripe phase I. In Fig[2£i we show the excitation spectrum 
in the phase II, calculated in the experimental conditions 



FIG. 2: (Color online) Lower branch of the excitation spec- 
trum in Phase II (a) (fi/^o = ^ g reen l me i ^/^o — 1-2, 
blue line) and in Phase III (b) (fi/fcjj = 2.25) as a function of 
the wave vector q. In phase II the spectrum is not symmetric 
and exhibits a roton minimum for negative q whose energy be- 
comes smaller and smaller as one approaches the transition to 
the stripe phase at 0,/k^ = 0.09. The dotted line corresponds 
to the usual Bogoliubov dispersion in the absence of spin-orbit 
and Raman coupling. The other parameters: Gi/kg — 0.12, 
7 = G2/G1 = 1Q- 3 . 



of [f|. In Fig. 2b we instead show the excitation spectrum 
in the phase III which exhibits symmetry by inversion of 
q into —q. The physical origin of the roton minimum is 
quite clear. In the phase II the ground state is degener- 
ate and it is very favorable for atoms to be transferred 
from the BEC state at p = fci to the empty state at 
p = —ki. The occurrence of the roton minimum is also 
reflected in a strong enhancement in the static response 
function x(q)- Notice that x(q), differently from u)(q), is 
always a symmetric function of q. The occurrence of the 
roton minimum in the excitation spectrum and the cor- 
responding enhancement of the static response represent 
a typical tendency of the system towards crystallization. 
We have investigated in detail the condition for the roton 
frequency being equal to zero (corresponding to a diver- 
gent behavior for x(<?))- A simple analytic expression 
for the corresponding value of the Raman coupling Q is 
obtained in the weak coupling limit G\,Gi <§C A)g were 
we find that the critical value exactly coincides with the 
value fli-u = 2k 2 ^/2"f / (1 + 27) characterizing the tran- 
sition between the plane wave and the stripe phases. For 
larger values of the coupling constants Gi and G2 we ex- 
pect that the critical value takes place for values of the 
Raman coupling smaller than the value at the transition, 
exhibiting the typical spinoidal behavior characterizing 
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the usual first-order liquid-crystal phase transition. 

In conclusion we have investigated the dynamic behav- 
ior of a Bose-Einstein condensate with spin-orbit cou- 
pling, pointing out the occurrence of novel features of 
high relevance for future experiments, like the strong 
quenching exhibited by the sound velocity near the 
second-order transition the plane wave and the single 
minimum phases, and the occurrence of a roton mini- 
mum in the excitation spectrum in the spin polarized 
plane wave phase. 

This work has been supported by ERC through the 
QGBE grant. 
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